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Metamaterial  Perfect  Absorber 
Analyzed  by  a  Meta-cavity 
Model  Consisting  of  Multilayer 
Metasurfaces 

Khagendra  Bhattarai1,  Sinhara  Silva1,  Kun  Song2,  Augustine  Urbas3,  Sang  Jun  Lee4, 
Zahyun  Ku3  &  Jiangfeng  Zhou  1 


We  demonstrate  that  the  metamaterial  perfect  absorber  behaves  as  a  meta-cavity  bounded  between  a 
resonant  metasurface  and  a  metallic  thin-film  reflector.  The  perfect  absorption  is  achieved  by  the  Fabry- 
Perot  cavity  resonance  via  multiple  reflections  between  the  "quasi-open"  boundary  of  resonator  and 
the  "close"  boundary  of  reflector.  The  characteristic  features  including  angle  independence,  ultra-thin 
thickness  and  strong  field  localization  can  be  well  explained  by  this  meta-cavity  model.  With  this  model, 
metamaterial  perfect  absorber  can  be  redefined  as  a  meta-cavity  exhibiting  high  Q-factor,  strong  field 
enhancement  and  extremely  high  photonic  density  of  states,  thereby  promising  novel  applications  for 
high  performance  sensor,  infrared  photodetector  and  cavity  quantum  electrodynamics  devices. 


The  metamaterial  perfect  absorber  (MPA),  exhibiting  nearly  unity  absorption  within  a  narrow  frequency  range1-6, 
has  attracted  tremendous  interest  recently  The  MPA  possesses  characteristic  features  of  angular-independence, 
high  Q-factor  and  strong  field  localization  that  have  inspired  a  wide  range  of  applications  including  electromag¬ 
netic  (EM)  wave  absorption3, 7> 8,  spectral5  and  spatial6  modulation  of  light,  selective  thermal  emission9,  thermal 
detection10,  infrared  photo  detection11-14,  circularly  polarized  light  detection15-17,  and  refractive  index  sensing  for 
gas18  and  liquid19-21  targets.  The  MPA  typically  consists  of  three  layers:  an  array  of  metallic  resonators  (e.g., 
cross-type  resonators9,22,  split-ring  resonators3  or  metallic  nanoparticles23  and  chiral  resonators15,24)  and  a  highly 
reflective  layer  (e.g.  metallic  film3, 9>  22, 23  or  metallic  mesh  grid3),  separated  by  a  subwavelength-thick  dielectric 
film  (spacer).  The  mechanism  of  MPA  has  been  explained  through  a  number  of  models  including  impedance 
matching3, 7’ 25  and  destructive  interference26.  In  the  impedance  matching  mechanism,  the  entire  three-layered 
structure  is  considered  as  a  thin  slab  made  from  a  homogeneous  medium  with  frequency- dependent  effective 
permittivity  E^uS)  and  effective  permeability  /iff(u;).  At  the  perfect  absorbing  frequency  (cj0),  the  effective  per¬ 
mittivity  and  permeability  reach  the  same  value  feeff(a;)  \U=UJ  =  /xff(u;)  \u=u ),  so  that  the  impedance  (both  real 
and  imaginary  parts)  matches  to  the  free  space,  ze{{(uj)  | u=  =  zef/(u0)  +  izeff"(uj0)  =  1.  Meanwhile,  the  effective 
refractive  index  «eff( uj)  at  u  =  exhibits  a  large  imaginary  part,  nef/f(u0)  %/(o;0),  thus  incident  EM  wave 
enters  MPA  without  any  reflection  and  then  rapidly  decays  to  zero  inside  the  MPA.  In  the  destructive  interference 
mechanism,  the  resonator  and  the  metallic  film  are  considered  as  two  decoupled  surfaces,  and  destructive  inter¬ 
ferences  of  waves  reflected  multiple  times  by  two  surfaces  lead  to  zero  reflection26, 27 

In  this  work,  we  use  a  transfer  matrix  method  to  assemble  the  MPA  from  three  functional  layers  and  obtain 
the  phase  and  the  amplitude  conditions  for  the  perfect  absorption.  Analyzing  these  conditions,  we  reinterpret 
the  MPA  structure  as  a  meta-cavity  bounded  between  a  “quasi-open”  boundary  of  a  resonant  metasurface  and  a 
“close”  boundary  of  a  metallic  film.  The  Fabry- Perot  (FP)  modes  are  achieved  when  the  phase  condition  is  sat¬ 
isfied,  which  leads  to  a  strong  resonance  with  non-perfect  absorption  (NPA).  However,  the  perfect  absorption 
(PA)  only  occurs  when  both  the  phase  and  the  amplitude  conditions  are  satisfied  simultaneously  Although  PA 


department  of  Physics,  University  of  South  Florida,  Tampa,  33620,  USA.  2Smart  Materials  Laboratory,  Department 
of  Applied  Physics,  Northwestern  Polytechnical  University,  Xi'an,  710129,  China.  3Materials  and  Manufacturing 
Directorate,  Air  Force  Research  Laboratory,  45433,  Wright-Patterson  AFB,  USA.  division  of  Convergence  Technology, 
Korea  Research  Institute  of  Standards  and  Science,  Daejeon,  34113,  Korea.  Khagendra  Bhattarai,  Sinhara  Silva  and 
Kun  Song  contributed  equally  to  this  work.  Correspondence  and  requests  for  materials  should  be  addressed  to  S.J.L. 
(email:  sjlee@kriss.re.kr)  orZ.K.  (email:  zahyun.ku.l.ctr@us.af.mil)  or  J.Z.  (email:  jiangfengz@usf.edu) 


SCIENTIFIC  REPORTS  |  7:  10569  |  D0l:10.1038/s41598-017-10520-w  1 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


www.nature.com/scientificreports/ 


Figure  1.  Effective  models  for  MPA.  Schematic  diagrams  of  (a)  metamaterial  perfect  absorber  (MPA)  and  its 
(d)  equivalent  single-layer  film  characterized  by  effective  <reff  and  /iff;  (b)  perfect  absorber  with  cross-wire, 
spacer  and  metal  ground  plate  considered  as  a  three-layer  structure  and  (e)  equivalent  three-layer  film  model 
where  the  cross-wire  is  considered  by  a  homogeneous  film  with  eeff  and  (i  ff .  (c)  View  of  perfect  absorber  along 
wave  propagation  direction.  The  MPA  is  surrounded  by  air  and  a  dielectric  substrate  with  optical  impedances  zsl 
and  zs2  at  two  interfaces,  respectively,  (f)  A  transmission  line  model  uses  a  two -port  network  S  to  describe  the 
MPA,  where  zY  and  z2  are  impedances  of  air  and  substrate,  and  zin  and  zout  are  input  and  output  impedances  of  S, 
respectively 


has  drawn  significant  attention,  NPA  are  more  commonly  used  in  practical  applications  such  as  sensors18-21,28 
and  detectors10-14.  Both  PA  and  NPA  can  be  well  explained  by  the  meta-cavity  model.  With  an  improved  retrieval 
method,  we  find  that  the  resonant  metasurface  operates  at  off-resonance  wavelengths,  thereby  acting  as  a  homog¬ 
enous  thin  film  with  high  dielectric  constant.  The  FP  model  redefines  MPA  as  a  meta-cavity  with  high  Q-factor 
and  extremely  high  photonic  density  of  states  (PDOS)  where  the  FP  modes  are  easily  tailored  by  the  geometric 
parameters  of  metasurface.  Our  work  can  pave  the  way  for  novel  photonic,  optoelectronic  and  cavity  quantum 
electro  dynamic  (QED)  applications29. 

Results 

Single-layer  effective  medium  model.  Our  perfect  absorber  consists  of  three  layers:  an  array  of  cross¬ 
wire  resonators,  a  dielectric  spacer  and  a  metal  ground  plane  (MGP)  as  shown  in  Fig.  la.  We  carry  out  a  3D  full- 
wave  simulation  to  solve  Maxwells  equations  and  obtain  numerical  solutions  by  Computer  Simulation 
Technology  Microwave  Studio  that  uses  a  finite  integration  technology30.  The  parameters  used  for  simulation  are 
as  follows:  period  of  cross-wire  array,  p  =  2  /ira,  length  and  width  of  the  cross-wire,  /  =1.7  /ira  and  w  =  0.4  /ira, 
permittivity  of  GaAs  substrate  and  a  dielectric  spacer,  eGaAs  =  11.56  and  es  =  2.28,  thickness  of  the  dielectric 
spacer,  ts  =  0.09  /ira,  and  thickness  of  cross- wire  and  MGP,  tc_w  =  tMGP  =  0.1  /ira.  The  gold  used  for  cross-wire 
and  MGP  is  described  by  the  Drude  model31  with  plasma  frequency,  ujp  =  1.37  x  1016  rad  Is  and  collision  fre¬ 
quency  ujc  =  4.08  x  1013  rad/s.  We  use  three  different  models  to  describe  the  perfect  absorber:  (i)  A  single-layer 
effective  medium  model  where  the  entire  structure  (Fig.  la)  is  considered  as  a  layer  of  homogeneous  medium 
(Fig.  Id)  characterized  by  the  effective  permittivity  (eeff)  and  effective  permeability  (/iff);  (ii)  A  three-layer  effec¬ 
tive  medium  model  where  the  cross-wire  resonator  (Fig.  lb)  is  considered  as  a  homogeneous  effective  film 
(Fig.  le)  with  eeff  and  /iff  on  top  of  two  layers  of  real  films,  the  spacer  and  MGP;  (iii)  A  transmission  line  model 
where  the  multiple-layer  optical  system  (Fig.  lc)  is  modeled  as  a  two-port  network  (Fig.  If)  with  input  impedance 
zin  and  output  impedance  zouV 

Since  the  MGP  prevents  the  EM  wave  from  propagating  through,  the  perfect  absorption  is  obtained  when  the 
reflection  reaches  zero  at  certain  wavelength  as  shown  in  Fig.  2d.  Among  three  models  shown  in  Fig.  1,  the 
single-layer  effective  medium  model  is  widely  used  in  various  MPA  works1,3,7,25,  where  the  electric  and  the  mag¬ 
netic  resonances  lead  to  the  impedance  matching  condition,  zeff  =  ^/iff/6gff  =  1  ,  at  the  perfect  absorption  wave¬ 
length.  Therefore,  the  incident  EM  wave  propagates  through  the  top  surface  of  effective  film  in  Fig.  Id  without  any 
reflection.  Meanwhile,  the  effective  film  exhibits  large  imaginary  effective  refractive  index,  n''ff  nerff,  so  that  the 
EM  wave  decays  exponentially  to  zero  before  it  exits  from  the  bottom  surface.  To  verify  this  interpretation,  we  use 
effective  medium  theory32  to  retrieve  the  effective  permittivity  ( eeff )  and  effective  permeability  (/iff)  from  the 
complex  transmission  and  reflection  coefficients,  t  and  r.  Due  to  the  asymmetric  structure  of  MPA  along  the 
propagation  direction  of  incident  light,  the  perfect  absorption  only  occurs  for  incident  illumination  from  the 
front-side  (cross-wire)  and  the  EM  wave  incident  from  the  back-side  (MGP)  will  be  totally  reflected.  Strictly 
speaking,  such  extremely  asymmetric  propagation  in  MPA  cannot  be  described  by  waves  traveling  in  a  simple 
medium  with  homogeneous  permittivity  and  permeability.  However,  since  we  only  consider  the  interaction  of 
MPA  with  incident  light  from  the  front-side  (cross-wire),  we  can  use  t  and  r  for  impinging  light  from  the  air-side 
to  obtain  eeff  and  /iff  as  shown  in  Fig.  2a, b,  respectively  We  find  a  Lorentzian  type  resonance  of  <reff  at  wavelength 
of  5.13  pm,  which  can  be  ascribed  to  the  first-order  dipole  resonance  mode  of  the  cross- wire  driven  by  the  electric 
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Figure  2.  Effective  parameters  for  single-layer  effective  medium  model.  Real  and  imaginary  parts  of  (a) 
effective  permittivity  and  (b)  effective  permeability  of  MPA  retrieved  by  using  a  single-layer  thin  film  model, 
(c)  Real  and  imaginary  parts  effective  impedance  calculated  by  single-layer  thin  film  model  (zeff,  solid  curves) 
and  effective  circuit  model  (zin,  symbols),  (d)  Simulated  reflection  spectra  of  the  actual  MPA  (blue-solid)  and 
structures  consisting  of  single-layer  (green- circle)  and  three-layer  (red-sold)  effective  films. 


field  of  the  incident  wave.  For  wavelength  A  >  6  /ira,  Re(eeff)  is  negative.  This  Drude  type  response  is  attributed 
to  the  MGP,  which  exhibits  the  typical  characteristic  of  permittivity  of  bulk  metal,  but  with  much  lower  plasma 
frequency  due  to  diluted  average  electron  density  resulting  from  low  volume-ratio  of  MGP  inside  the  MPA  unit 
cell33.  The  effective  permittivity  /iff  also  shows  a  Lorentzian  type  resonance  at  wavelength  of  5.94  pm  driven  by  the 
magnetic  field.  This  magnetic  resonance  is  induced  by  the  anti-parallel  currents  on  the  cross-wire  and  the  MGP, 
and  can  be  qualitatively  explained  by  the  method  of  image  theory7.  The  solid  lines  in  Fig.  2c  are  the  real  (red  color) 
and  the  imaginary  (blue  color)  parts  of  the  effective  impedance  zeff,  which  show  strong  dispersive  response  as  a 
function  of  wavelength. 

At  the  perfect  absorption  wavelength  X  =  5.94  pm  (feff  =  2.432  +  24.175/  and  /iff  =  1.003  +  25.517/),  the 
effective  impedance  is  found  to  be  zeff  =  1.02  +  0.03/,  which  matches  with  the  impedance  of  air  (zair  =  1)  very 
well.  The  effective  refractive  index,  neff  =  ^ce ff  •  /iff  =  1.736  +  24.849/  has  a  large  imaginary  part  (i.e., 
n"  =  24.849).  These  calculated  effective  constitutive  parameters  very  well  support  the  single-layer  effective 
medium  model.  Specifically,  the  matched  impedance  allows  the  incident  wave  passing  through  the  top  surface  of 
the  effective  thin  film  (Fig.  Id)  without  any  reflection.  Meanwhile,  the  intensity  of  the  wave  is  I  =  I0e~2az ,  where 
a  =  nk  =  —  26.28  \im~l  (i.e.,  incident  wave  decays  at  a  rate  of  He  for  every  0.019  fim  distance  traveling 

inside  the  mAv).  When  the  incident  wave  is  reflected  by  MGP  and  travels  a  round-trip  inside  MPA  back  to  the 
front- surface,  its  intensity  reduces  to  I/I0  =  e~Aa'ts  =  7.8  x  10-5.  Thus  99.998%  of  the  wave  energy  is  absorbed 
inside  the  MPA.  As  shown  in  Fig.  2d,  the  second  absorption  peak  can  be  also  found  at  X  =  1.79  pm  and  the  effec¬ 
tive  impedance  zeff  =  1.046  +  0.003/  is  found  to  be  very  close  to  air  (zair  =  1)  again.  Other  effective  parameters 
take  the  following  values:  eeff  =  0.3008  +  7.157/,  jli{{  =  0.378  +  7.8331/ and  neff  =  0.338  +  7.49/.  To  the  best  of 
our  knowledge,  this  second  absorption  peak  (X  =  1.79  pm)  has  not  been  thoroughly  studied  in  the  past.  As  will  be 
shown  later,  it  results  from  the  second-order  dipole  resonance  of  the  cross-wire.  The  effective  impedance  of  the 
MPA  can  also  be  obtained  by  the  effective  circuit  model34.  Figure  If  shows  an  effective  circuit  for  the  MPA,  where 
the  MPA  is  described  as  a  two-port  element  S  in  connection  with  loads  zx  (air)  and  z2  (substrate).  We  calculated 
the  input  impedance  zin  based  on  transmission  line  theory  and  plotted  it  as  triangle  (real  part)  and  circle  (imagi¬ 
nary  part)  symbols  in  Fig.  2c.  Although  zin  was  calculated  with  complete  different  methodology  (see  details  in 
Methods  section),  it  shows  excellent  agreement  with  zeff  obtained  from  single-layer  thin  film  model.  Using  the 
retrieved  effective  permittivity  and  permeability  eeff(cj)  and  /iff(u;),  we  performed  full-wave  simulation  for  the 
single-layer  effective  thin  film  (Fig.  Id).  Figure  2d  shows  that  the  simulated  reflection  (green-circle)  completely 
overlaps  the  reflection  from  the  original  absorber  structure  (blue-solid  line)  as  expected.  The  single  layer  effective 
medium  model  can  provide  fundamental  insight  for  MPA  (e.g.,  reflection  and  transmission),  however  this  model 
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Figure  3.  Effective  parameters  for  three-layer  effective  medium  model  and  perfect  absorption  conditions. 

Real  (red)  and  Imaginary  (blue)  parts  of  (a)  effective  permittivity  (edi)  and  (b)  effective  permeability  (/iff).  (c) 
|r12|  (red)  and  \a  •  r23|  (blue):  terms  used  in  the  amplitude  condition  (Eq.  2).  Inset  depicts  various  transmission 
and  reflection  coefficients  of  MPA  structure,  (d)  0  (red),  7  (blue),  2(3  (green),  </>(r21)  (black)  and  0(r23)  (orange): 
terms  used  in  the  phase  condition  (Eqs  3, 4).  Two  perfect  absorption  regions  at  5.94  pm  and  1.79  pm  are 
highlighted  in  green. 


cannot  provide  the  microscopic  mechanism  of  the  absorption,  cannot  very  well  explain  NPA  and  can  miss  some 
important  features  such  as  ultra-thin  thickness. 

Three-layer  effective  medium  model.  To  better  understand  the  underlying  mechanism  that  the  wave 
travels  and  decays  inside  the  MPA,  we  have  studied  the  transmission  and  reflection  of  waves  at  each  constituting 
layers  and  used  a  transfer  matrix  method  to  obtain  the  overall  reflection/absorption  properties  of  MPA.  We  con¬ 
sidered  the  cross-wire  as  a  metasurface  and  described  it  as  a  homogeneous  thin  film  with  effective  permittivity  eeff 
and  effective  permeability  /iff.  The  effective  parameters  were  calculated  using  transmission  and  reflection  coeffi¬ 
cients  of  an  air/cross- wire/ spacer  configuration  through  an  improved  retrieval  method  to  take  account  of  the 
asymmetric  structure  due  to  different  incoming  and  outgoing  medium  (air  and  spacer35).  As  shown  in  Fig.  3a,  the 
effective  permittivity  of  cross-wire  shows  a  Lorentzian  type  electric  resonance  at  wavelength  of  Xr  =  4.75  /ira, 
which  is  shorter  than  the  electric  resonance  obtained  by  the  single-layer  model  (5.13  /ira).  The  resonance  at 
4.75  /ira  is  purely  electric  since  the  effective  permeability  (Fig.  3b)  takes  a  constant  value  of  Re(/iff)  =  0.8  except  a 
typical  anti-resonance  at  Xr  due  to  the  periodicity  effect36.  At  wavelength  of  1.59  /ira,  we  found  the  second-order 
electric  dipole  resonance  mode  of  the  cross-wire,  which  is  much  weaker  and  has  negligible  influence  on  the  value 
of  eeff.  With  those  values  of  eeff  and  peff,  we  used  a  thin  film  of  the  same  thickness  to  replace  the  cross-wire  in  the 
absorber  structure.  3D  Full- wave  simulation  of  the  three-layer  thin  film  model  (Fig.  le)  was  carried  out.  The 
simulated  reflection  (red-solid  curve)  is  plotted  in  Fig.  2d,  which  matches  very  well  with  simulation  of  actual 
MPA  structure  (blue- solid)  at  the  first  absorption  wavelength  (A  =  5.94  /ira)  and  is  slightly  off  at  the  second 
absorption  wavelength  (A  =  1.79  /ira).  We  ascribe  the  small  discrepancy  at  A  =  1.79  /ira  to  possible  effect  of 
MGP  on  the  resonance  of  cross-wire,  which  makes  slightly  noticeable  difference  of  reflections  for  wavelengths 
below  2.5  /ira.  At  two  absorption  peaks,  eeff  reads  value  of  64.72  +  7.08i  and  4.19  +  1.28/,  and  /iff  is 
0.78  —  0.018/  and  0.8  —  0.018/,  respectively.  The  role  of  the  cross-wire  can  be  completely  understood  as  a  thin 
film  with  these  eeff  and  /iff.  In  spite  of  minor  discrepancy  at  short  wavelengths,  the  excellent  agreement  between 
the  three-layer  thin  film  model  and  actual  MPA  proves  that  the  cross-wire  layer  and  the  MGP  layer  are  decoupled, 
i.e.  neither  the  existence  of  MGP  will  affect  the  resonance  of  cross- wire  nor  the  resonant  current  on  the  cross-wire 
will  affect  MGP  where  the  latter  only  plays  a  role  of  a  reflecting  mirror.  Therefore,  the  transfer  matrix  method  as 
a  common  approach  to  study  multiple-layered  structure  can  be  used  to  obtain  the  overall  property  of  MPA  from 
each  layer  and  to  investigate  the  absorption  mechanism  of  EM  waves  inside  the  MPA.  The  overall  reflection  of  the 
three-layer  structure  can  be  calculated  by  multiplying  the  transfer  matrix  of  each  layer,  M  =  Ml  •  M2  •  M3,  as 
given  below  (more  details  can  be  found  in  our  previous  works27, 35). 
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r  _  rl2  +  ar23e 

1  —  r2lr23e  ^  (1) 

In  Eq.  1,  ru  and  r21  are  the  reflection  coefficients  of  the  cross-wire  from  front  (air)-  and  back  (spacer) -side,  respec¬ 
tively.  r23  is  the  reflection  coefficient  of  the  MGP.  (3  =  ns  •  k  ■  ts  is  the  propagating  phase  in  the  dielectric  layer,  where 
ns,  ts  and  k  are  the  refractive  index,  the  spacer  thickness  and  the  wave  vector  in  free  space,  respectively. 
a  =  hJn  ~  r2iri2>  where  *12  and  *21  are  the  transmission  coefficients  through  the  cross-wire  along  forward  (air / 
cross- wire/spacer)  and  backward  (spacer/ cross- wire/air)  directions,  respectively  Although  a  is  strictly  equal  to  1  at 
the  interface  between  two  homogenous  media,  a  ^  1  for  the  cross-wire  around  the  resonance  wavelength.  This  is 
caused  by  the  structural  asymmetry  depending  on  the  propagation  direction  of  incident  light.  For  perfect  absorption 
(i.e.,  r  =  0)  according  to  Eq.  1,  we  obtain  the  following  conditions  for  the  amplitude  and  the  phase,  respectively: 

Vu\  =  l«  '  r23\  (2) 

9  =  4>(r12)  -  </>(r23)  -  <j>(a )  -  2(3  =  (2 n  +  l)7r,  \n\  =  0,  1,  2,  ....  (3) 

The  coefficients  r12,  r21,  f21  and  tl2  are  obtained  from  simulation  of  the  air/cross-wire/spacer  configuration  and  r23 
is  obtained  from  the  simulation  of  the  spacer/MGP  configuration. 

Around  the  perfect  absorption  wavelength  \  =  5.94  pm,  both  the  amplitude  condition  in  Eq.  2 
(|r12|  =  \a  •  r23|)  and  the  phase  condition  ( 6  =  7r)  are  satisfied  simultaneously  as  shown  in  Fig.  3c  and  d,  respec¬ 
tively  For  wavelengths  below  5.94  pm,  9  decreases  as  the  wavelength  decreases  and  finally  reaches  — ty  at 
X  =  1.79  pm.  Meanwhile,  the  amplitude  condition  is  also  satisfied  at  X  =  1.79  pm,  thereby  leading  to  the  second 
absorption  peak.  By  carefully  examining  Fig.  3d,  we  find  that  the  amplitude  condition  (|r12|  =  | a  ■  r23|  at 
X  =  5.61  pm)  and  the  phase  condition  (9  =  tt  at  X  =  5.88  pm)  are  not  satisfied  at  exact  same  wavelength.  This 
explains  the  fact  that  the  reflection  shown  in  Fig.  2d  does  not  reach  zero  completely 
(min(  | r|2 )  =  Rmin  =  3.68  x  10-4  at  X  =  5.94  pm).  The  phase  condition  in  Eq.  3  suggests  a  Fabry-Perot  cavity- like 
mechanism.  The  cross-wire  acts  as  a  “mirror”  that  reflects  light  due  to  its  high  effective  permittivity  and  forms  a 
cavity  together  with  the  MGP.  The  incident  wave  enters  the  MPA  and  is  reflected  multiple  times  by  MGP  and 

cross-wire.  The  FP  cavity  resonances  establish  when  the  round-trip  phase  is  fulfilled  with  the  following 

condition: 


7  =  </>(r21)  +  </>(r23)  +  2 (3  =  2nm ,  \m\  =  0,  1,  2,  ....  (4) 

where  </>(r21)  and  </>(r23)  are  the  phases  of  reflection  coefficients  at  the  interfaces  of  the  spacer/cross-wire  and  the 
spacer/MGP  (Fig.  3d),  respectively  When  7  equals  to  2mn,  the  waves  that  travel  multiple  round-trips  between 
cross-wire  and  MGP  are  in  phase  and  hence  interfere  constructively,  thereby  leading  to  resonant  cavity  modes.  As 
the  wavelength  decreases,  7  (blue  curve  in  Fig.  3d)  increases  monotonically  and  reaches  the  values  of  0  and  2n  at 
wavelength  of  5.94  pm  and  1.82  pm,  respectively  (matching  very  well  with  two  perfect  absorptions  wavelengths). 
As  shown  in  Fig.  3d,  the  main  contribution  to  7  comes  from  </>(r21),  which  increases  from  0.897T  to  2.637T  as  the 
wavelength  decreases  from  5.94  pm  to  1.82  pm.  The  significant  change  of  </>(r21)  is  caused  by  6eff,  which  changes 
from  64.72  +  7.08i  to  4. 19  +  1.28/  in  the  same  wavelength  range  as  shown  in  Fig.  3a.  Note  that  the  phase  condi¬ 
tion  of  Eq.  4  only  guarantees  the  FP  cavity  modes,  resulting  in  the  absorption  peaks.  To  achieve  perfect  absorp¬ 
tion,  the  amplitude  condition  in  Eq.  2  needs  to  be  satisfied  at  the  same  wavelength  as  well.  Nevertheless,  Eqs  3  and  4 
well  explain  the  NPA  peaks  of  a  variety  of  resonator/ spacer/MGP  structures  that  have  been  widely  used  in  sensing 
and  detecting  applications10-14, 18-21,28  where  the  absorption  peaks  do  not  necessarily  reach  100%  (i.e.,  R  =  0).  The 
ultra-thin  thickness  of  MPA  can  be  also  explained  by  the  cavity  model.  As  shown  in  Fig.  3d,  at  the  first  perfect 
absorption  wavelength  X  =  5.94  pm  where  7=0,  the  phase  </>(r23)  «  —  ty  and  </>(r21)  is  slightly  less  than  7r,  i.e., 
</>(r21)  =  7r  —  <5  with  6  <C  ir.  To  achieve  7  =  0,  it  is  required  [3  «  <5/2,  thus  the  spacer  thickness  ts  is  found  to  be 
6\/4n  <  A.  In  our  simulations,  we  obtained  6  =  0. 1  l7r  and  thus  ts  =  0.0275  •  As  «  0. 108  fim,  which  is  close  to 
the  actual  thickness  ts  =  0.09  fim  «  A5/44  where  As  =  A lns  =  3.94  fim  being  the  wavelength  inside  the  spacer 
medium.  This  thickness  is  extremely  thin  in  comparison  to  a  regular  FP  cavity  with  size  of  \SI2.  Note  that  the  key 
fact  that  leads  to  this  ultra-small  size  of  cavity  is  </>(r21)  being  slighly  smaller  than  n  due  to  high  effective  permittiv¬ 
ity  of  the  cross-wire  layer. 

Another  important  characteristic  feature  of  MPA  is  the  angular  independence  of  the  incident  wave.  This  is 
very  useful  in  thermal  imaging/sensing  and  energy  harvesting  applications  to  receive  off-normal  incident  lights. 
We  performed  simulations  of  actual  MPA  structure,  single-layer  and  three-layer  thin  film  model  with  the  polar 
angle  9  (in  spherical  coordinates  as  displayed  in  the  inset  of  Fig.  4)  of  incident  wave  varying  from  0°  to  80°.  The 
reflection  spectra  are  shown  in  Fig.  4a-c  for  TE  wave  and  Fig.  4a-f  for  TM  wave,  respectively  As  shown  in  Fig.  4a, 
the  wavelength  of  absorption  peak  at  X  =  5.94  pm  for  TE  wave  does  not  shift  as  the  incident  angle  increases  from 
0°  to  80°,  showing  good  angular  independence.  Reflections  of  both  the  single-layer  (Fig.  4b)  and  the  three-layer 
(Fig.  4c)  models  well  match  with  that  of  actual  MPA  structure  (Fig.  4a).  For  TM  wave,  the  reflection  of  the  MPA 

is  blue-shifted  as  9  is  increased  by  Aa  =  X80°  ~~  x°°  ~  4%  as  shown  in  Fig.  4d.  Both  the  single-layer  and  the 

|  ^0°  | 

three-layer  models  also  show  the  blue-shift  with  Aa=7.4%  and  9.1%,  respectively,  which  is  slightly  higher  than 
actual  MPA.  The  three-layer  model  can  qualitatively  explain  the  angular  dependence  for  TE  and  TM  modes.  The 
cross-wire  layer  exhibits  an  electric  resonance  excited  by  the  electric  field  of  incident  wave.  For  TE  mode,  the 
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Figure  4.  Color  plot  of  angular  dependence  ( 6 )  of  reflection.  Reflection  for  TE-polarized  incidence  for  (a) 
actual  MPA  structure,  (b)  single-layer  thin  film  model  and  (c)  three-layer  thin  film  model,  (d)-(f)  Reflection  for 
TM-polarized  incidence  for  corresponding  structures. 


electric  field  (See  the  inset  of  Fig.  4a)  is  always  parallel  to  the  cross-wire  in  spite  of  changing  the  polar  angle  9. 
Thus,  the  resonance  of  cross-wire  does  not  change  with  incident  angle.  This  leads  to  the  angular  independence  of 
absorption  for  the  entire  MPA  structure.  For  TM  mode,  however,  the  electric  field  direction  changes  with  0  (see 
the  schematic  diagram  as  shown  in  the  inset  of  Fig.  4d),  resulting  in  shifting  the  wavelength  of  absorption  peak. 
At  short  wavelengths  (A  <  p),  incident  waves  become  significantly  diffracted  through  the  cross- wires  and  lead  to 
shifting  and  splitting  of  the  second  absorption  peak  at  \=1.79  pm  as  indicated  in  Fig.  4a  and  d.  These  diffractions 
do  not  exist  for  two  effective  film  models  as  shown  in  Fig.  4b,c,e  and  f.  In  addition,  we  observed  new  absorption 
peaks  at  around  2.75  pm  and  4.75  pm  (Fig.  4d)  for  TM  wave,  which  becomes  more  significant  for  6  >  50°.  These 
absorption  peaks  (reflection  dips)  result  from  resonance  modes  of  cross-wires  and  MGP  induced  by  the  vertical 
component  of  electric  field  (Ez),  which  is  nonzero  for  TM  wave  and  is  increased  with  9. 

Higher  order  FP  modes.  Equations  3  and  4  indicate  that  the  phase  conditions  are  satisfied  when  (3  increases 
by  an  integer  multiple  of  27t,  thereby  leading  to  a  series  of  FP  cavity  modes.  The  increment  of  2n  in  (3  corresponds 
to  an  increase  in  the  thickness  of  spacer  (ts)  by  Xs/2.  For  the  first  absorption  peak  at  X  =  5.94  pm  shown  in  Fig.  2d, 
the  increment  of  thickness  is  calculated  as  As/2  =  1.97  fim  and  the  corresponding  n- th  order  thickness  is  given 
by  tsn  =  tsl  +  n{Xs/2)  where  tsl  =  0.09  fim  is  the  thickness  of  the  original  MPA.  Figure  5a  shows  the  reflection 
spectra  for  the  three  lowest  order  thickness,  tsl  =  0.09  fim,ts2  =  2.06  yumandfs3  =  4.03  fim.  The  reflection  spec¬ 
tra  for  the  thickness  ts2  =  2.06  fim  (red  curve)  and  ts3  =  4.03  fim  (blue  curve)  show  the  absorption  peaks  at 
wavelengths  of  5.65  pm  and  5.81  pm,  respectively,  which  are  well  matched  with  the  absorption  peak  at  wavelength 
of  5.94  pm  for  the  original  MPA  (i.e.,  tsl  =  0.09  fim  for  1st  order  thickness).  We  attribute  the  minor  discrepancy  in 
wavelengths  to  the  slight  changes  of  cross-wire  resonances  due  to  different  separation  between  cross-wire  and 
MGP.  Figure  5b  shows  the  reflections  of  three-layer  thin  film  model  using  the  same  thickness.  The  results  very 
well  match  the  reflections  of  actual  MPA  (Fig.  5a).  We  recognize  these  resonances  as  the  three  lowest  order  FP 
cavity  modes,  which  is  confirmed  by  observing  the  x-component  of  electric  field  Ex  inside  the  MPA  structure.  The 
right-side  figures  in  Fig.  5c-e  show  the  real  part  of  Ex  along  z-direction  inside  the  spacer  at  the  wavelength  of 
absorption  peaks,  where  Ex  is  evaluated  along  a  straight  line  in  z-direction  below  the  center  of  cross-wire  (from 
z  =  0,  the  bottom  of  cross-wire  to  z  =  tsn,  the  spacer  thickness).  For  the  first-order  thickness  tsl  =  0.09  fim, 
Re(iy  decreases  monotonically  from  positive  value  to  zero  as  z  increases  from  0  (at  the  cross-wire/spacer  inter¬ 
face)  to  0.09  fim  (at  the  spacer/MGP  interface).  Such  Ex  distribution  results  from  the  boundary  conditions.  At  the 
spacer/MGP  interface,  the  high  conductivity  of  MGP  forms  a  “close”  boundary  where  the  tangential  electric  field 
is  forced  to  be  zero  (EUx  =  0).  At  the  cross-wire/spacer  interface,  the  high  effective  permittivity 
(ee ff  =  64.72  •  (1  +  0. 11/)  at  \=5.94  pm)  of  cross-wire  provides  a  “quasi-open”  boundary  which  allows  non-zero 
tangential  electric  field  (Et  x  ^  0).  As  the  thickness  of  spacer  increases  to  ts2  =  2.06  fim,  a  half-period  of  sine 
curve  is  added  to  Ex  field  for  the  first  order  thickness  (Ex  as  shown  on  the  right  side  of  Fig.  5c)  as  shown  in  Fig.  5d, 
which  results  in  forming  the  second-order  FP  cavity  mode.  Similarly,  Ex  for  ts3  =  4.03  fim  (Fig.  5e)  shows  the 
fundamental  mode  plus  a  full-period  of  sine  curve,  leading  to  the  third-order  FP  cavity  mode.  The  three-layer 
thin  film  model  gives  the  same  Ex  profiles  as  shown  in  Fig.  5f-h.  All  the  Ex  fields  for  actual  MPA  (Fig.  5c-e)  and 
for  effective  thin  film  model  (Fig.  5f-h)  exhibit  constant  phases  (not  shown  here)  along  z-direction,  thereby  indi¬ 
cating  standing  waves  as  we  expected  according  to  the  abovementioned  FP  cavity  analysis.  The  left-side  figures  in 
Fig.  5c-e  show  the  electric  current  distributions  on  the  cross-wire  and  MGP  for  the  three  lowest  orders  of  FP 
mode.  In  all  three  figures,  the  currents  on  the  cross-wires  show  a  half-wavelength  of  cosine  curve  along  the 


SCIENTIFIC  REPORTS  1 7:  10569  |  D0l:10.1038/s41598-017-10520-w  g  6 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


www.nature.com/scientificreports/ 


Figure  5.  Reflection,  Ex  field  and  current  distributions  for  lowest  three  orders  FP  modes,  (a)  Simulated 
reflection  spectra  for  actual  MPA  structures  with  the  spacer  thickness  of  tsl  =  0.09  fim  (1st  order:  black), 
ts2  =  ts i  +  As/2  =  2.06  fim  (2nd  order:  red)  andf53  =  tsl  +  As  =  4.03  fim  (3rd  order:  blue),  (b)  The 
corresponding  reflection  spectra  for  three-layer  thin  film  model  with  spacer  thickness  tsl  =  0.09  /ira, 
ts2  =  2.06  fim  and  ts3  =  4.03  fim.  Top  views  of  current  density  on  the  cross-wire  and  MPG  for  (c)  1st  order, 

(d)  2nd  order  and  (e)  3rd  order  MPA.  Insets  placed  on  the  right  side  of  (c)-(e)  show  the  ^-component  of  electric 
field  (Ex)  along  z-direction  (propagation  direction).  Current  density  on  the  meta-film  representing  cross-wire 
and  the  MGP  for  (f)  1st  order,  (g)  2nd  order  and  (h)  3rd  order.  Insets  on  the  right  side  of  (f)-(h)  display^  along 
z-direction  for  the  three-layer  thin  film  structures. 


horizontal  direction  with  maximum  at  the  center  and  zero  at  ends.  This  reveals  a  first-order  dipole  resonance 
mode  of  the  cross-wires  and  matches  with  the  Lorentzian  resonance  of  se{{  shown  in  Fig.  3a.  Figure  5c  shows  the 
current  on  the  MGP  flowing  toward  the  opposite  direction  with  the  maximum  current  in  the  center  region.  The 
anti-parallel  currents  on  the  cross-wire  and  MGP  in  Fig.  5c  induce  a  strong  magnetic  moment  that  is  responsible 
for  the  magnetic  resonance  for  single-layer  thin  film  model  shown  in  Fig.  2b.  We  also  observe  the  same 
anti-parallel  currents  in  the  simulation  for  the  three-layer  thin  film  model  as  shown  in  Fig.  5f.  Such  current  dis¬ 
tributions  can  be  well  explained  by  the  Ex  profile  shown  in  Fig.  5f.  The  tangential  electric  field  Et  x  induces  the 
electric  current,  }x  =  iu;eEt  x ,  where  the  e  for  cross-wire  and  MGP  are  given  by  eeff  and  <fMGP  (using  the  Drude 
model31)  respectively  At  wavelength  of  X  =  5.94  pm,  the  real  parts  of  e  for  cross-wire  and  MGP  have  the  opposite 
signs,  i.e.,  Re(^)  >  0  and  Re(<rMGP)  <  0.  For  the  first-order  cavity  shown  in  Fig.  5c, f,  Ex  in  metafilm  and  MGP 
are  in  phase,  and  thereby  result  in  anti-parallel  currents.  For  the  higher  order  cavities  (e.g.  Figure  5g,h  for  the 
second-  and  third-order  cavities,  respectively),  when  ts  is  increased  by  As/2,  the  phase  of  Ex  in  MGP  increases  by 
7r  and  the  current  in  MPG  changes  to  the  opposite  direction. 

We  expect  the  FP  cavity  model  for  MPA  to  lead  to  novel  applications.  As  one  example,  the  MPA  provides 
much  stronger  local  electromagnetic  fields  enhancement  than  typical  plasmonic  resonators  due  to  its  large 
Q-factor.  Such  local  field  enhancement  is  always  demanded  to  improve  the  performance  of  sensor/detector  at 
infrared  and  THz  regimes11, 21, 37 .  However,  the  region  between  resonator  and  MGP  in  the  MPA,  where  the  fields 
are  mostly  enhanced,  are  too  thin  to  fit  most  sensing  and  detector  devices.  With  the  cavity  model,  we  can  design 
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the  higher  order  cavities  that  can  fit  in  the  required  thickness  of  practical  devices.  For  instance,  infrared  focal 
plane  array  device  working  at  6  pm  wavelength37,  is  composed  of  an  active  layer  (20  stacks  of  InGaAs/AlGaAs 
quantum  wells)  and  top/bottom  contact  layers  (n-doped  GaAs).  The  total  thickness  is  about  2  pm  which  is  much 
larger  than  the  first-order  MPA  thickness  of  0.09  pm ,  but  is  very  close  to  the  second-order  thickness  of  2.06  pm. 
Furthermore,  this  thickness  is  tunable  by  varying  the  size  of  the  resonator  (cross-wire)  to  fit  the  requirement  of 
existing  devices.  Another  potential  application  is  to  use  the  high  photonic  density  of  states  (PDOS)  for  cavity 
quantum  electrodynamic  (QED)  applications29, 38, 39.  PDOS  is  inversely  proportional  to  the  volume  of  a  cavity  The 
first-order  cavity  of  MPA  has  thickness  of ~  A/44,  much  smaller  than  typical  FP  cavity  (A/2).  This  size  is  also  much 
smaller  than  typical  plasmonic  resonators29  which  EM  fields  normally  extend  to  distance  of  a  few  wavelengths. 
Our  calculation  shows  the  MPA  cavity  can  improve  the  PDOS  by  ~80  times  as  compared  to  plasmonic  resonators. 
Such  high  PDOS  promises  great  potential  in  cavity  QED  applications29, 38, 39. 

Discussion 

In  summary,  we  developed  a  multi-layered  metasurface  model  for  MPA.  With  a  transfer  matrix  analysis,  we 
found  a  FP  cavity  mechanism  that  models  the  MPA  as  a  meta-cavity  bounded  between  the  “qusi-open”  bound¬ 
ary  of  resonator  and  the  “close”  boundary  of  MGP.  The  FP  cavity  model  well  explains  the  characteristic  features 
including  ultra-thin  thickness  and  angular  independence.  We  also  found  the  higher  order  FP  cavity  modes  when 
the  thickness  of  MPA  (spacer)  increases  by  multiples  of  a  half- wavelength.  The  strong  field  enhancement  of  cavity 
resonances  and  high  PDOS  promise  novel  applications  in  sensor/detector  and  cavity  QED  device. 


Methods 

The  effective  impedance  zinof  the  circuit  model  is  calculated  from  the  reflection  coefficients  Tin  and  Tout  at  the  the 
input  and  output  termial  of  S,  using  the  following  equations34: 


p  _  c  |  ^12^ouAl 

in_  11  +  Y^7s~ 

i  r2°22 


(5) 


^out  —  ^2: 


^12^X1^21 

1  -  qSn 


(6) 


where  Sn,  S12,  S21  and  S22  are  the  S-parameters,  which  can  be  obtained  from  simulations  of  entire  MPA  structure. 
The  coupled  equations  (Eqs  5  and  6)  can  be  solved  analytically: 
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We  then  solved  the  input  and  output  impedances  from  Tin  and  Tout: 

=  1  +  rin 

in  1  -  rin  (9) 


1  +  r0i 

1  -  D, 


(10) 


For  MPA,  we  only  consider  the  input  impedance  zin  as  shown  in  Fig.  If  for  the  incident  wave  from  the  front 
side. 
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